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Computing  Machinery  Meeting  at 
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Arnold  o.  Mengel 

*4 

An  analogue  computer  of  the  KLAC  type,  although! 
designed  primarily  for  simulation,  is  a  valuable  tool  in 
apnlied  research.  Jeverul  modifications  have  increased 
the  flexibility  and  usefulness  of  the  HAhl  HLAC . 


The  FiEAC  proved  its  worth  in  the  study  of  the 
application  of  the  calculus  of  variations  to  the  optimiza¬ 
tion  of  aircraft  flight  paths.  Prior  to  the  installation 
of  the  HEAC  several  reports  on  the  theoretical  aspects  of 
the  problem  were  written,  but  reached  no  conclusions  because 
of  the  complexity  of  the  equations.  HEAC  solutions  were  of 
a  sumrising  nature  that  not  only  led  to  un^-xoectedly 
practical  answers,  but  also  suggested  revisions  of  the 
theory. 
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The  calculus  of  variations  na-e  beer  s  r^-rH  c d  ajfc=tfefctrD  to 


CX  /•» 


answer  'Such  questions  as 

^ — { 1  )  hnt  path  should  a  missile  with  a  fixed  period  of  thrust 
fly  to  maximise  its  range  at  a  specified  final  velocity 
and  altitude® 


(^(^)  shat  path  should  an  aircraft  fly  after  take-off  to 
minimize  its  time  of  flight  to  level-flight,  combat 
velocity  at  a  specified  altitude** 


The  resulting  equations  are  computationally  difficult  to 
handle  because  of  their  complexity  and  quasi -stable  nature.  More¬ 
over,  the  new  variables  introduced  by  the  calculus  of  variations 
have  no  apparent  Dhysical  meaning,  v/hich  makes  the  analysis  of  their 
influence  difficult. 

(j  o  f~&  ‘-vv 

This  paper  outlines  the  experiences^ vw  'hrwe  'nad-at^ia^-tr 
"solving" “'tal'cnr itrs'*0'f 'va rt  ATT o ns  pro bl  em  s  on  the  REAC.  A  modified 
form  of  the  equations  (presented  as  an  appendix  to  this  paper)  was 
developed  which  not  only  was  more  satisfactory  computationally,  but 
also  showed  that  the  calculus  of  variations  enuations  described  the 
motion  of  a  body  similar  to  the  one  under  study.  The  nature  of  the 
IlEAC  solutions  suggested  that  the  steady-state,  or  mid-path,  trajec¬ 
tory  and  the  transient  trajectories  from  the  steady-state  path  to  the 
end-points  could  be  computed  separately.  The  results  using  this  method 
agreed  well  with  the  ■*irexactM>  solution^  and  gave  a  tre  i  w.  A  v. a  savings 
in  computing  tine*  sir^ce  much  of  the  trio  1-and-error  process  of 
meeti-ng  end-conditionv*>««4iJie  usual  process  was  eliminated  and  the 
stability  problem  was  'n.  mPr  r. »  tm- ^  trajectories 


-•‘Quotation  marks  have  been  placed  around  the  word  exact  ^trecause 
certain  approximations  were  necessary  in  the  computation. 
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Consider  a  body  of 
weight  w  fly  n g  _n  two-dim¬ 
ensional  space  with  a  velocity 
v  in  a  direction  making  an 
arrle  0  with  the  x-axis  as 
shown  in  Figure  1.  The  total 

forces5'*  acting  on  the  body  4 

are  resolved  into  four  com-  / 

patients,  the  weight  vector  w,  / 

7w  /  At  V 

a  force  due  to  the  change  of  X  / 

weight  -  V ,  a  f^rce  vector 

M  acting;  along;  the  axis  of  the  ^ 

body,  and  a  force  K  normal  to  r  1^** 

its  axis,  with  ff  =  7?  +  Tf.  ^ - - - 

Th  >  angle  of  attack  <*■  is  the 

angle  between  the  velocity  Figure  1.  Force-  on  the  Body 

vector  V  and  the  axis  of  the 

body.  The  rnarnitudes  of  M  and  N  are  assumed  ti  b«  functions  c 


the  altitude  y,  velocity  v,  and  angle  of  attack  ot  .  Defining- 
If  as  the  position  vector  c-f  t he  hody,  with 

R  =  /x~<  + 

the  dynamic  equations  are 
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with  initial  conditions  * 

X  -  X,  ,  y  -  y,  ,  V  »  V,  ,  9  -  9,  ,  w  -  w,  ,  and  t  -  t,. 

The  problem  is  to  determine  the  rath  with  terminal  conditions 

X  -  X2f  y  "  y^»  v  "  »  9  “  9fcl  w  »  w2,  at  t  «  t2 

such  that  one  of  the  terminal  variables  is  a  minimum  or  maximum. 
The  problem  may  be  modified  by  not  specifying  certain  terminal 
values,  but  if  v2  or  ©^  are  not  prescribed  (or  optimized)  the 
problem  becomes  singular.  The  problem  may  be  modified  further 
by  specifying  constraints. 


DISCUSSION 

The  problem  posed  in  the  previous  section  is  not  treated 
in  the  standard  texts  on  the  calculus  of  variations1 ,  and  the  first 
treatment  of  the  problem  at  RAND  was  by  E.  M.  Liebhold  .  Later, 

Dr.  Magnus  Hestenes,  as  a  consultant,  developed  a  more  complete 
and  elegant  presentation  2*3#.  The  work  of  this  paper  was  done 
with  the  cooperation  of  Dr.  Hestenes,  Kenneth  Martin,  and  Rogeh^Snow. 


Application  of  the  calculus  of  variations  to  the  above 
problem  yields  a  set  of  Euler-Lagrange  differential  equations  (and 
variables)  to  be  solved  simultaneously  with  the  dynamic  equations 
of  the  system  under  study.  The  initial  conditions  of  the  Euler 
variables  may  be  varied  to  cause  the  trajectory  to  pass  through  any 
possible  set  of  end-conditions.  Certain  conditions  must  be  satis¬ 
fied  to  assure  the  solution  is  an  optimum^. 
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Unfortunately  the  Euler  equations  are  usually  considerably 
more  complicated  than  the  dynamic  equations.  Furthermore,  the 
adjustment  of  the  initial  values  of  the  Euler  variables  to  meet 
specified  end-conditions  is  a  tedious  t rial-and-error  process  made 
more  difficult  by  the  absence  of  any  physical  significance  for  these 
variables.  To  make  matters  worse,  the  equations  turn  out  to  have 
a  quasi-stable  nature  very  similar  to  that  of  a  neophyte  tight-rope 
walker.  This  combination  of  woes  led  to  the  development  of  the 
equations  appearing  in  the  appendix. 

The  remainder  of  this  paper  will  be  a  blow-by-blow 
description  of  our  trials  and  tribulations  with  the  calculus  of 
variations.  We  hope  to  continue  work  on  the  problem  at  some  future 
time. 

Problem  1 

Our  first  experience  with  a  calculus  of  variations  problem 
came  when  our  acquaintanceship  with  the  REAC  was  just  starting  and 
nearly  resulted  in  a  parting  of  the  ways. 

The  purpose  of  this  problem  was  to  test  how  sensitive  a 
maximum  was  to  approximations  mad*  id  .the  Euler  equations.  The 

body  studied  had  constant  weight,  noliltr*  constant  thrust,  *£4 
drag  proportional  to  the  square  of  the  velocity.  The  resulting 
dynamic  and  Euler  equations  were  quite  simple. 

Our  faces  were  crimson  when  the  solutions  of  the  approx¬ 
imate  Euler  Equations  gave  maximums  several  percent  greater  than 
those  given  by  the  exact  equations.  No  manner  of  kicking,  balancing, 
swearing,  checking,  and  even  cheating  a  bit  on  settings  could  deviate 
the  REAC  from  its  opinion  that  the  approximate  solution  was  the 
better. 


P-199 
1-15-51- 
Pare  5 


Finally  in  desperation  I  derived  the  Euler  equations  and 
was  overjoyed  to  find  the  problem  had  been  submitted  with  a  wrong 
sign  in  one  of  the  Euler  equations.  Thus  the  KEAC  not ' only  was 
vindicated,  but  at  the  same  time  was  made  a  hero. 

The  corrected  equations  denonstrar  eu ,  a.  .honed,  that 
the  optimums  are  flat  enough  that  approximate  trajectories  give 
results  for  all  practical  purposes  as  good  as  the  optimum  trajectories. 

Problem  2 

This  Droblem  was  a  sequal  to  the  above  -  an  attempt  to 
get  a  "feel"  for  trajectory  optimization  and  to  find  the  influence 
of  the  initial  values  of  the  Euler  variables.  A  plot  v:as  made  of 
final  velocity  v2  as  a  function  of  final  range  x2  for  the  body  of 
Problem  1  at  specified  final  values  of  altitude  yA  .  and  time  t^  with 
zero  final  vertical  velocity  (©2  *  0  or  180°). 

It  was  hoped  that  a  linearized  set  of  equations  would 
help  in  establishing  the  relationships  between  the  initial  values 
of  the  Euler  variables  to  give  the  proper  end-conditions,  but  since 
the  partials  changed  rapidly  and  the  inaccuracies  introduced  by 
taking  differences  were  large,  a  plot  of  two  of  the  variables 
versus  the  arctan  of  the  third  proved  more  helpful  and  at  the  same 
time  gave  some  insight  as  to  their  influence  on  the  system.  The 
form  of  the  plot  of  *2  versus  x2  is 
illustrated  in  Figure  2.  Point  A  is 
the  terminal  pair  for  maximum  x2  with 
x2  not  specified,  while  Point  8  is 
the  terminal  pair  for  maximum  x2  with 
x2  not  specified.  Actually,  only  the 
arc  length  designated  AS  is  of  any 
practical  interest.  It  is  apparent 

.that  since  the  lower  portion  of  the  Figure  2.  Values  of 

x^  versus  x2  for 
specified  y2  ,  ©2,  and  t2 
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curve  gives  minimum  values  of  while  the  upper  nortion  gives 
the  maximum  values,  all  other  oossible  terminal  conditions  lie  within 
the  closed  curve.  For  example,  the  system  can  not  have  a  terminal 
condition  such  as  specified  by  Point  P. 

Problem  3 

To  maximize  the  range  of  a  missile  it  is  necessary  to 
study  the  powei — on  and  nower-off  (or  glide)  stages  separately. 

This  problem  was  the  power-on  trajectory  optimization.  /.  set  of 
values  similar  to  the  arc  AS  of  Figure  2  was  desired  for  several 
values  of  92  and  yk  and  was  to  be  used  as  the  initial  values  of 
the  power-off  trajectories. 

The  Euler  eouations  for  the  power-on  trajectory  were  too 
complicated  to  fit  on  our  RE AC  and  a  t rial-and-error  process  was 
required  to  give  pseudo-optimum  paths.  Since  the  peaks  appeared 
very  flat  and  no  better  hand-controlled  paths  could  be  found,  use 
of  the  resulting  trajectories  seemed  justified. 

The  flight  path  was  first  programmed  by  letting 

• 

9  =  a  +  Tj-t.  The  value  of  a  was  varied,  b  and  9]  adjusted  to  meet 
the  end-conditions  and  the  sets  found  that  gave  maxi mums.  The 
angle  of  attack  giving  the  programmed  9  was  computed  and  used  in 
evaluating  thrust  and  drag.  A  second  method  of  nrogramming  with 

9  =  0  for  0  <  t  <  T 

0  =  a  +  Tj(t  -  T)  for  1  <  t  <  t„ 

gave  nearly  identical  maximums  with  what  appeared  at  the  time  as 
more  realistic  trajectories. 
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Problem 


During  the  power-off  portion  of  the  trajectory  the 
velocity  remained  well  above  the  spepd  of  sound,  the  coefficients 
of  lift  and  drag  were  reasonably  well  behaved  and,  as  a  consequence, 
it  was  possible  to  compute  the  dynamic  plus  Euler  equations  on  the 
KEAC .  Since  it  was  desired  to  maximize  range  with  no  regard  to 
time,  letting  range  be  the  independent  variable  and,  eliminating 
time  greatly  simplified  the  computation. 

The  trajectories  were  very  sensitive  to  changes  in  the 
initial  conditions  of  the  two  Euler  variables  and  many  trial  nuns 
had  to  be  made  to  locate  the  region  of  interest.  A  typical  set 
of  trajectories  wi th  the  initial  value  of  one  of  the  Euler  variables 
fixed  and  the  other  varied  is  shown  in  Figure  3*  A  charge  of  the  initial 
value  by  about  5  per  cent  would  swing  the  trajectory  from  oath 
A  to  B.  This  great  sensitivity  and  the  quasi-stable  condition  is 
understandable  when  one  considers  the  many  end  conditions  that  must 
be  possible  no  matter  how  undesirable  they  may  be  (i.e.,  vertical 
flight,  reversed  headings,  etc.) 
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It  ;.'u s  aonan-r.t  that  the  path  giving  maximum  range  v;a c 
between  and  3.  Paths  C  and  )  v:er>-  next  found  and  differed  by 
the  order  of  about  one  oer  cent  in  the  initial  value  of  the  dale r 
variable.  The  optimum  path  was  between  t:  em;  so  using  initial 
values  between  theirs  we  found  the  paths  £,  anti  ?,  having  the  order 
of  0.1  per  cent  difference  in  their  initial  conditions.  This 
discussion  is  based  upon  the  assumption  that  the  initial  value 
held  constant  was  at  the  proper  value  for  the  optimum  path. 
Actually,  the  initial  valua  of  the  sec  ond  Lult-r  variable  would  have 
to  be  "jockeyed"  along  with  the  other  to  keep  the  solutions  from 
"blowing  up".  3y  the  time  path  E  and  F  were  reached,  the  cirift  in 
the  d.c.  amplifiers  prevented  a  continuation  of  the  above  procedure 
since  the  changes  in  the  initial  conditions  became  of  the  same 
order  of  magnitude  as  the  drift.  f  simple  change  of  scale-factor 
v/ as  impossible  since  the  variables  range  1  from  +100  to  -100  volts 
despite  their  sensitivity  to  initial  conditions. 

An  initial  value  bounded  by  those  of  naths  E  and  F  was 
selected,  the  problem  stopped  at  the  node  P  and  readings  made  of 
the  variables.  These  values,  with  perturbations  on  the  sulor 
variables,  were  used  as  new  initial  conditions.  This  process  was 
continued  until  the  final  velocity  dropned  to  vk  and  require  i  a 
full  day  of  computing.  Oince  this  process  had  to  be  reoeated  for 
several  sets  of  v,,  9,,  and  Vj ,  the  prospects  looked  grim. 

However,  inspection  of  the  first  two  days’  results  le  : 
to  a  much  simpler  computation.  First  it  was  noticed  that  the 
oaths  oscillated  about  th«  altitude  for  which  lift  equalled  weight 
at  the  angle  of  attack  C<c  maximizing  lift  ever  drag.  Moreover,  the 
optimum  paths  approach#?  i  this  critical  altitude  ir>  much  td 
way  an  underdamre  i  servo-mechanism  res;  Uids  to  a  step  function. 
Hence,  we  tried  orogramrni ng  (X  as 

&  ■  (X  ( 1  -  K  sin  o] 
c 
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where  Q(  was  a  function  of  velocity  and  K  was  varied  until  a 
c 

maxi  man  range  was  obtained.  Tb.e  unstable  appearing  end-condit  ions 
could  be  obtained  by  changing  the  sign  in  front  of  K;  by  varying 
K  and  the  time  of  reversal  of  sign  any  end-conditions  could  be 
realized  by  paths  similar  to  A  through  F  of  Figure  3»  The  values 
of  range  obtained  by  this  approximate  method  were  within  0.2  per 
cent  of  those  found  by  the  "exact"  solution,  and  required  about  one 
minute  of  REAC  time  compared  to  the  eight  hours  required  for  the 
"exact"  solution. 

Problem  5 

The  Euler  and  dynamic  equations  for  an  aircraft  flying 
below  the  speed  of  sound  are  simple  enough  to  handle  on  our  REAC. 

In  this  problem  we  attempted  to  minimize  time  of  flight  from 
take-off  to  level  flight,  combat  velocity  at  a  specified  altitude. 

The  first  attempt  at  solution  showed  that  the  equations 
were  even  less  stable  than  those  of  Problem  L, .  In  fact,  the  paths 
dove  into  the  ground  no  matter  how  the  Euler  variables  were 
adjusted. 

At  this  time  the  work  of  the  appendix  was  started  and  it 
became  apparent  that  the  constraint  y  >  0  must  be  added,  since 
an  initial  diving  path  to  pick-up  velocity  was  optimum.  Obtaining 
a  solution  even  with  the  computationally  improved  equations  of  the 
exact  form  in  the  appendix  was  still  more  difficult  than  for  those 
of  Problem  U  and  again  we  were  forced  to  develop  an  approximate 
method  . 

Previous  work  indicated  that  if  the  initial  and  final 
values  were  sufficiently  separated,  the  mid-paths  of  all  trajectories 
were  nearly  the  same  except  near  the  transitions  to  the  boundary 
values.  Hence,  as  shown  in  the  appendix,  a  quasi-steady  state 
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solution  v; as  found  giving  the  optimum  midpath  for  all  cases.  The 
comDlete  set  of  equations  was  used  only  for  the  transition  from 
the  mid-path  to  the  end-points.  Also  equations  were  found  giving 
either  the  ontiraum  velocity  or  rate  of  climb  at  which  to  fly  as 
functions  of  altitude,  and  the  computation  of  the  steady-state 
trajectory  became  simple  enough  for  hand  solution. 

The  dynamic  equations  for  the  hypothetical  aircraft 
were 

x  *  v  cos  ©  ft/sec 
y  ='v  sin  0  ft/ sec 

v  *  [l6  -  1.36»10“5  v2  -  3  •  1  3  *1 0-3  i  va. )  2  J<T  -  32.2  sin  9  ft/ sec2 

©  *  3.l3*10-3crva  -  5  ^  radians/sec. 

■  relative  air  density  «  e  a^. 

Figure  4  shows  the  form  of  the  quasi-steady  state 

solution  plus  one  initial  set  of  conditions  and  several  final 

conditions.  If  the  final  velocity  is  approximately  ten  per  cent 

or  more  than  the  climbing  velocity,  a  final  dive  (as  in  Path  A)  is 

found  to  give  quite  a  saving  in  time  over  a  oath  leveling  off  at 

the  desired  altitude  and  flying  level  until  the  soecified  velocity 

is  reached.  A  final  steep  climb  before  leveling-off  is  proper  when 

a  velocity  less  than  the  critical  climbing  velocity  is  desire!.  The 

initial  portion  of  the  trajectory  was  found  by  flying  level  until 

•  • 
ana  was  called  for  making  ©  >  0  at  which  time  the  constraint  6  >  0 

was  removed,  as  shown  in  Figure  5»  A  more  satisfactory  method  would 

involve  coming  off  the  steady-state  oath  with  the  comnlet^  equations 

as  in  the  final  portion  of  the  trajectories  but  with  * i me  running 

backwards  ? 

Figure  6  and  the  equivalent  nomogram  of  Figur°  7  give  not. 
only  the  minimum  time  of  flight  for  specified  values  of  v^ ,  , 

and  ©^ ,  but  also  the  maximum  value  of  with  v k and  specified 
and  the  maximum  value  of  v^  with  t2,  y^ ,  and  specifier. 


studied 


< 


where  cr 
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The  optimum  velocity  and  rate  of  climb  for  the  hypothet¬ 
ical  aircraft  to  fly  us  functions  of  altitude  are  given  in  Figures 
$  and  9  respectively.' 

Comparisons  of  or, her  flight  programs  with  the  optimum  are 
illustrated  by  Figures  10  and  11.  The  aircraft  climbed  at  a  velocity 
Av  greater  than  the  optimum  in  the  first  case  and  at  a  constant 
velocity  in  the  second  case.  In  both  cases  when  the  climbing 
velocity  was  less  than  the  specified  final  value,  the  aircraft 
leveled  off  at  the  specified  altitude  and  flew  level  until  the 
desired  velocity  was  reached.  For  those  few'  cases  where  the  climbing 
velocity  was  greater  than  the  final  value,  the  aircraft  pulled  into 
a  climb  before  leveling-off. 

It  is  difficult  to  justify  mathematically  the  use  of  the 
quasi-steady  state  path  for  the  major  portion  of  the  flight  oath 
and  the  use  of  the  exact  baler  equations  only  for  the  transition 
from  the  steady  state  oath  to  and  fr  n  the  end  points.  However, 
the  fact  that  the  exact  soluti  in:  follow  the  steady  state  path 
smoothly  for  a  while  and  then  oscillate  about  the  steady  state 
path  when  they  "blow-up”  is  some  justification.  If  then  all 
optimum  paths  have  very  nearly  identical  midpaths,  it  becomes 
clear  why  the  quasi-stable  nature  of  th.i  equations  is  necessary  to 
make  it  possible  to  meet  .all  end-roints* 
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APPENDIX 


BASIC  EQUATIONS  (assuming  ^  »  C) 

De  f iningv 

H  =X*7  >  fl-V* 

where  A  and  p.  are  variable  multipliers,  the  Euler-Lagrange  equat 
optimizing  the  trajectory  described  by  Figure  1  and  Equations  1 

-  -  VrK  1 

P  *  -  VVH  \  (2) 

Ha  =  C 


where 


V, 


rradient  of  H  in  K  soace 


R 

— V 

=  gradient  of  H  in  V  space 


(the  symbol^  may  be  viewed  as  a  partial  gradient,  i.e.,  the 
gradient  when  all  variables  but  R  are  held  fixed}.  Notice  that 
the  dynamic  equations  of  motion  are 


R 


H1 

(1 


V  -VhH 


(la) 


Letting  the  argument  of  ^  be  $,  define 
p  -  (J>  -  9  -  a 
G  * 


M  cos  p  +  N  sin  p  -  •Fr  J 


(3) 


or  G  is  the  projection  of  on  the  vector  "jT.  Then, 


*The  advantages  of  the  vectorial  derivation  were  pointed  out  By 
Roger  Snow.  T;  is  appendix  assumes  the  reader  is  familiar  with 
the  theory  of  the  calculus  of  variations  given  in  RI->1C'Q» 


ions 

are 


{ 


H  =  X"v  +  £  (UG  ♦  M 


A  -  -  £  f*V  G 

w  d 

i7  *  -■  S  ^vvG  -  A 

H  -  |£  -  0 

C(  W'  (Jot 


>  (4) 


Treating  p  as  an  independent  variable 
Gv  =  Mv  cos  p  ♦  Nv  sin  p  ' 


x 


0 


G  =  M  cos  p  +  N  sin  P 
y  v  r  y  r 

Ga  =  Ma  cos  p  +  sin  p 
Gp  ■  -  M  sin  p  +  N  cos  p 


f  151 


Since  p  *  (J>  -  0  -  a 


G-  =  -  n, 


'0 


JP 


and 


dG  r  r 

aa  Ga  GP 


H*  *  ^(Ga  -  0p)f-0 
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yields 

Ga  "  Gp 

As  a  consequence,  the  last  two  equations  of  set  5  yield 
Ma  cos  p  +  Na  sin  p  =  -  M  sin  p  +  N  cos  p 
or 


N  -  Ma 

e  *  Cretan  --- 


(6) 


P-199 
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Reducing  the  vector  equations  1  and  4  to  polar  form  gives 
x  =  v  cos  9 

y  »  v  sin  9 

•  v  =  cos  ot  -  N  sin  a  )  -  g  sin  © 

v©  =  sin  a.  +  N  cos  a  )  -  g  cos  © 

P  =  ^[-Gv  cos  (a  +  P)  ♦  qp  sin  (<*+P)]p-  A  sin  $  -  Ax  cos 


pi  -  g[« 


Gv  sin 


(a+p)  +  -2  cos  (cx+p)Jpi  -  A  cos  $  Ax  sin  (J) 


A  =  -  £  G  p 
'v  w  v r 


N  -  M, 


Cretan 


in  which  G^  has  been  replaced  by  the  enuivalent  -  0C 

It  can  be  seen 

that  the  Euler  equations 

describe  a  system  acted 

upon  by  the  forces  shown  ^ 

in  Figure  12.  Thus  these  \  /\s<V 

equations  can  be  considered  ^  \  & 

as  describing,  the  dynamics  * 

of  an  "Euler — craft"  having  m 

an  axial  drag  force  uG„,  N' 

G  v  V  V 

a  normal  lift  force  |u_»  ,  ^  f 

and  two  pseudo-gravity 
forces  Av  (variable)  and 
K  (constant).  The  axis 

A 

of  the  "Luler-cra ft"  lies 

along  the  V  vector  and  its  Figure  12.  Forces 

-velocity"  vector  TT  is  "Kuler-craft" 

(«♦  P)  radians  above  that  of  its  axis. 


Figure  120  Forces  on  the 
"  Eu  1  er-cra  ft" 
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» P  PROX I  A  T  L  L '  ■  U  A  TUNS 


In  p-eneral  the  forces  or  an  aircraft  are  piven  by 
S  K  ■  T(y,v)  -  f i  ly,v)<rv‘  1 


w 


w 


N  =  f  ^  ( y  ,  v )  cr  v  *  a 


(8 


wher  e 

<r  ~  e  is  the  relative  air  density  and  f,  and  f2  are 
nearly  constaiit  for  velocities  belov:  the  speed  of  sound.  Ac  a 
consequence 


p  =  arctan 


K  -  Ma 
!77  7  W 


irctan- 


f^  cr  v^  a 


T  —  f  j  cr  v^  +  f  ^  o' v^ 
since  T  -  ^CTv2  «  f^  CTv^  and  O.  is  small. 

If  we  neglect  the  influence  of  the  axial  forces  on  the 

•  * 

equations  for  0  and  0,  and  let  p  =  ,  sin  c*  =  3.,  sin  2C*-  =  2°t, 

cos  a  -  1 ,  and  cos  20 (  =  i ,  equations  7  become 

x  =  v  cos  0 

y  =  v  sin  0 

v  =  ^ ( M  -  fi<*)  -  g  sin  0 
w 


arctan  a  ~  a 


© 


£  N  g  cos  0 


w  v 


M  _  g 

/u  w 


[ - C M v  ♦  <*Nv)  +  (Ma  ♦  -  A  sin  0  -  Ax  cos  0  ^  (9) 

id~ 


0  =  §  [wa  +  aNa]  *  cos  0  +  sin  ^ 


A  =o 

X 
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Since  Q  =  l©  +  2a),  o(N  =  N ,  vh 

’  »  ’  v 


-N  (  f ,  «  constant )  , 


and  Ma  =  0,  the  Holer  equations  can  be  modified  to 

2a  =  ^  COS  © - ^  cos  (9  +  2a)  +  -^2  sin  (0  +  2Ct) 

jr  =  -  £  Kv  "  £in  (&  +  2 a]  ~  -j r  cos  +  2*)  >  do) 

-A*  =  -  +  at;  1 

m  w  L  y  yJ 


If  the  final  range  is  rot  specified,  A  —  0 ,  £>  =  — 

X  n 

is  a  convenient  change  of  variable,  and  the  Euler  equations  become 


2<x  = 


cos  0 


-  E,cos  (©  +  2a' 


a  =  -  £f?v;  +  cxn 

^  Wl  y 


•y]  *  it.  *.:«)♦§ 


/ 

Notice  that  once  v,  ,  0,,  and  a,  are  set,  the  setting  of  ieter- 
mines  «, ,  thus  clarifying  the  influence  of  the  Euler  variables. 

QUASI-STEADY  oVATr  CuNDITliNS 

During  the  quasi-stcadv  ntat  :  condition,  v,  0,  a,  a,  and 
^will  be  small  and  we  can  use  the  anoroxi motions 


g  CO 


S  0  ar  £  N 
w 


g  sin  e 


u;  -  Na; 


£,  c;  H 

V 


f  Fm  +  <XN  "1  sin  (0  +  2 CH '  +  E  ;;  1 

wl  y  yJ  V  Lv  w  vj 


Since 


_n  9  «  ^[M  *  Not] 


g  sin  (0  ♦  .M]  5;  (g  cos  9)20.  + 


'T 
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or 


ast  ecus  t 

i  o  n 

bee  .mes  ^ 

'id.  & 

_ 

♦  oci,')  =  (M  -  ai) 

.  g  a>y 

c?  v 

v  /I — 

a 

_ dL)  q 

=  F  -"v  = 

H  jd  t3 

of  change  of  energy 

(  £ 

<9y 

a  v 

1  ■ 

equut ion 

so  lv 

ed  for  ot 

s imultaneous 

1  y  v:  i  t  h  t ;  .r  dynamic 

it.  i  ons 


yields  the  optimum  quasi-irteady  state  trajectory. 

Putting  equation  13  in  the  form  of  set  d  and  assuming 
that  the  derivatives  of  the  f  functions  with  resoe^t  to  y  and  v 
have  only  a  second-order  effect,  v:°  find 

—  (t  ♦  af,  crv2  -  af ^(va)^  -  T  ♦  f ,crv*  -  f*cr(va)*  *  vTy  -  2f',crv^ 

a  y  v 


or 


*  ty  -  vTv  -T  *  (3  *  iziiilli 


(U) 


£  C 1  +  &~-)  cos  9 


After  take-off  the  plane  should  fly  level  (with  9  >  0  constraint) 

•  • 

until  9  becomes  positive  lat  which  tine  9  >  .0  constraint  is  removed). 

While  this  equation  is  suitable  for  HE AC  computation,  in 
practice  either  the  optimum  velocity  or  optimum  rate  of  climb  as 
function-s  of  altitude  would  be  sore  usef.l  to  a  pilot.  Moreover, 
the(rate  of  climb)  versus  altitude  curve  can  be  integrated  with 
respect  to  altitude  tc  find  the  tine  of  rlight  alone  the  quasi¬ 
steady  state  path.  The  variable  a  can  be  eliminated  by  approximating 
it  by  that  value  it  should  have  to  keep  9=0.  To  find  the  optimum 
velocity  or  rate  of  climb  v:e  must  solve  the  simultaneous,  equations 

E2(1  ♦*—)  cos*1  9*  f^v-p^.  Ty  -  vTv  -  T  +  (3  +a^)f,  ^  v^] 
oiUd  _  S  ^  f,  _  /T  -  f,  cry^Y 

g  TprT?  L  \  e  )  J  * 


sin  9  ~ 


T  -  f. 


(15) 
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As  a  second  joproach ,  and  £,  may  be  eliminated  by 
reducing  the  two  first  order  differential  equations  of  set  11  to 
one  second  order  equation  of  the  form 

&  ♦  i  v y , v , o ) <>  ♦  r’t l y , v , r? }0(  =  :'jiy,v,b  .  (lo) 

Durinp  the  quasi-steady  state  portion  of  the  fiipht  path  tnis 
equation  reduces  to 


0(  = 

"  FY(y,v,S'/  * 

is  a  third  approach,  we  c-: 


(17) 


p  cos  b 


-  £,003  (b 


to  solve  implicitly 6  for 


w  L1%V  +  WNyJ  *  ^  sin  +  A’ 


a)  +  -  1’.  | 
v;  vj 


to  solve  implicitly  for  a(, 

£(M  -  CXI,  )  -  p  sin  b  *  0 
to  solve  implicitly  for  b,  and 
“  N  -  p  cos  b  =0 

to  solve  for  implicitly  fcr  v.  This  method  works  particularly  v/ell 
on  the  itEAC  for  it  is  possible  to  wire  the  nroblem  such  that  one 
can  switch  from  the  steady-state  solution  to  the  transient  solution 
durinp  a  run.  It  also  is  convenient  to  allow  reversal  of  time  to 
find  initial  as  well  as  final  transient  naths  ? 
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